Stability of bound states in the light-front Yukawa model 



M. Mangin-Brinet, J. Carbonell 
Institut des Sciences Nucleaires, 53, Av. des Martyrs, 38026 Grenoble, France 

V.A. Karmanov 
Lebedev Physical Institute, Lemnsky Pr. 53, 117924 Moscow, Russia 
(Dated: February 1, 2008) 

We show that in the system of two fermions interacting by scalar exchange, the solutions for 
jT_Q+ bound states are stable without any cutoff regularization for coupling constant below some 
critical value. 
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I. INTRODUCTION 

A promising approach to solve the bound-state problem in field theory is the Light-Front Dynamics (LFD) ^, ||, |j. 
In this approach, the state vector is defined on the surface t + z — 0. The bound states in the Yukawa model (two 
fermions interacting by scalar exchange) were studied in papers ^ using Tamm-Dancoff method. It was found in 
particular that because the dominating kernel at large momenta tends to a constant, the binding energy of the 7 = 0"*" 
state is cutoff dependent, what requires the renormalization of Hamiltonian. 

The two-fermion wave functions were also considered in the explicitly covariant version of LFD 0. The state vector 
is there defined on the plane given by the invariant equation uj-x — with uj'^ — 0. The wave functions for the deuteron 
1^ and the pn scattering = 0+ state were found in a perturbative way and successfully applied to calculating 
the deuteron e.m. form factors |10| recently measured at TJNAF In papers ||^, |^ the NN interaction was taken 
from the Bonn model which includes scalar, pseudoscalar and vector boson exchanges with cutoffs in the A^A'^-meson 
vertices. The cutoff dependence of the spectrum was not analysed. LFD equations have now been solved exactly for 
a two fermion system in the ladder approximation with different boson exchange couplings [ p^ . 

In the present paper, we consider the solutions of the LFD equations for the scalar exchange only. We investigate 
the stability of the bound state solutions relative to the cutoff, disregarding the self energy contribution and renormal- 
ization. We have found a critical phenomenon for the cutoff dependence of the energy levels. The J = 0+ solutions 
are stable - not cutoff dependent when it tends to infinity - for coupling constant below some critical value. On the 
contrary, for coupling constants exceeding the critical value, the system is unstable (collapses) even in absence of the 
most singular interaction kernel in the Hamiltonian matrix. Our results are compared to those obtained in pi. 



In section II we write down the equations in a convenient form for analysis. In section [II, the cutoff dependence 
of the energy levels is studied analytically. In section the results of numerical calculations are given. Section ^ 
contains the concluding remarks. 



II. EQUATIONS 

We start with the system of equations for the two spin components $1 and $2 which contain a superposition of 
singlet and triplet spin states: 

{M^-M^)<i>iiR^,x) = £2 J [VniR±,x;R'^,x')^i{R'^,x') + Vi2{R±,x;R'^,x')<p2{R'^,x')]R'^dR'^dx', 

{M^ -Ml)<^2{R±,x) = j [V2i{RA_,x-R'^,x')'^i{R'^,x') + V22{R±,x-R'^,x')<^2{R'±,x')]R'^dR'^^^^ (1) 

with a = 1^, Mq — and kernels Vij given by: 

V^j = / v^j[{K'^ + i?)x{l-x)x'{l~x')]-^d(t)' 
Jo 

Vll = R±R'j^{2xx' ~x~ x') + [R]_x'{1 - x') + R'lx{l ~x)-m^(x + x'){2 ~x~ x')] cos0', 
V12 = m [Ra_{x + 3a;' - 2xx' - 2x''^) - R'^{x' + ix- 2x'x - 2x^) cos^'] , 
V21 = m [i?l(x' + 3x- 2x'x -2x^)- R±{x + 3x' - 2xx' - 2x'^) cos0'] , 
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V22 = R±R'j^{2xx' - X - x')cos(f>' + [R^^x'il- x') + R'lx{l - x) - m'^{x + x'){2~ X - x')] . (2) 
The momentum reads (see eq. (7.27) from [Q): 



i^2=„^2:^^_- + _ 2R^R'^ COS 0' + ^R'l + (x^ - x) ^ ' ^^1, - (3) 

X \ T I X X' \X'[\ — X') J 

for X < x' and with the replacements a; <-> x' R± ^ i?^ for a; > a;'. Functions $i_2 sue normahzed as: 

J {\<i>iiR^,x)\^ + \^2iR±,x)\''} R^dR^dx ^ I. (4) 

Equations (||) and kernels (||) are taken from |^ (eqs. 3.1a-b and C1-C4) with the notations = = $2 

and A; = R±, q = R'^, y = x' . 

In view of further analysis it is useful to introduce the new functions /i_2 given by: 

with N{x)='^-^^x{l-x] ~ ^ /-i?, m 



and the variables k^O {k > Q, Q < 9 < tt) defined by: 

i?^ = fcsin0, .= i('l--^^V (5) 
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In this new representation, the normalization condition (0) obtains the form 



d^k 

{fl{k,e) + fl{k,e)] — = i, (6) 



{2nf J eu 



where £k = ym^+k^, and one gets for fi the system of equations: 

f d^k' 

[^{e + m')-M^]h{k,e) = -— / [Ki^{k,e-k\e')h{k',e') + K^2{k,e-k',e')h{k\e')] 



[A{e + m^)-M^]h{Ke) = J [K2i{k,0;k',e')h{k',0') + K22{k,0;k',0')h{k',e')]—. (7) 

For convenience, we keep in (|7|) the integration over dcj)' , though it has been already performed in Kij. The kernels 
Kij are linearly expressed in terms of Vij and their analytical expressions read: 



27r 



(if 2 + ^2)^2£^£^, 2n ' 

Kii = -Q!7r [2fc2fc'2 + Sfc^m^ + 3fc'W + 4to'^ - 2kk'ekek' cos6lcos6l' - A:/c'(fc2 + fc'^ + 2™^) sin 6* sin 0' cos 0'] , 
K12 = — Q!7rm(A:^ — fc'^) (fc' sin 0' + fc sin0 cos (/>') , 

K21 = -Q!7rm(fc'^ - fc^) (fcsin6' + fc'sin6''coS(/>') , (8) 
K22 = -aTT [(2fc2fc'2 + 3fc2m^ + 3A;'^m2 + Am'^ - 2kk'ek£k' cos6'cos6'') cost/)' - fc/c'(fc2 + fc'^ + 2m^) sin0sin6''] , 

and (see eq. (3.60) from 0): 



if^ = fc2 + fc'2 - 2fcA:' f 1 + — cos e cos 6*' - 2kk' sin 6* sin 6*' cos </>' + f e?. + e?, ~-M'^ 



k cos A;' cos 0' 



£k 



In the variables fc, 0, the kinetic energy in (|^) is quadratic on k, the kernels are smooth in 0, and the stability of the 
binding energy is related to the kernels behavior at large k. 
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In the explicitly covariant version of LFD, the states are labeled by the eigenvalues J corresponding to the appropri- 
ate angular momentum operator Q J = Jq + Jfi, containing, besides the free operator Jp, the term J^j = —i[n x dn], 
where in cm. -system n = w/|w|. The two fermion wave function with J = is determined by two spin components 
The J = 1 wave function is determined by six components j^] and the equations are split in two subsystems, 
distinguished by the eigenvalues a = 0,1 oi A = (J-n)^. The J = 1, a = subsystem includes two components and the 
J = 1, a = 1 one includes four. The coupled equations obtained in this way (two for J = and two-|-four for J = 1) 
correspond to the 2-1-2-1-4 systems from In this classification, the equations (|^) with the kernels are written 
for the J = state, corresponding to (l-t-,2— ) from In the explicitly covariant LFD these equations are directly 
obtained in the form (^. Their direct derivation, together with the J = 1 case, will be given in a more detailed 
publication [ |T2[ . The relation between $i and fi and the corresponding relation between Kij and Vij result from the 
different representations of the spinors used in g and 0|. Below we present also the results for the J = l,a — state, 
whose components are related by a linear combination to <i)i^,$^+ defined in ||]. To distinguish the states, when 
necessary, we will attach to the kernels the index J = or J = 1 and omit a = 0. 

We would like to emphasize that the equations ([T]) solved in Q are related to our explicitly covariant LFD equations 
only by a linear tranformation of the components and a variable change. Both equations are thus strictly equivalent. 



III. THE CUTOFF DEPENDENCE OF THE BINDING ENERGY 



We consider the equations on the finite interval < k < k„iax- The dependence of the solution on the cutoff 
kmax in the limit fcmaa: — > oo is determined by the kernels asymptotics. Let us first analyze the kernel Kn. In the 
(fc, fc')-plane, when both k,k' oo with a fixed ratio fc'/fc = 7, this kernel tends to a constant. From the expressions 
we find the asymptotics: 

2 r ^Ani0,0'n), if7< 1 



Kn^-^{ An{e,e\l/j) .^^^^ (9) 



TO 



with 



, ^ 27(1 - cos 61 cos 6>0 - (1 + 7=^) sin 9 sin 9' cos (j) 

' ~ ^ Jo 2^ (1 + 72) (1 I cos 61 - cos9'\ - cos 9 cos 9') - 2-f sin 9 sin 0' cos (f> ' ^ ' 

and a' = a/{2mn). For convenience, we extracted the factor ^/7 in eq. (^). When k' is fixed and k — > 00, Kn 
decreases like 1/fc, and analogously for k fixed, fc' — > 00. 

We will compare the kernels in eq. (Q) with the kernel corresponding to the non-relativistic potential U{r) = —a'/r^. 
To get the exact correspondence, we take the S-wave Schrodinger equation for the wave function il){k), with the 
kernel U(k,k'), Fourier transform of U{r) = —a'/r'^, and we substitute there ip{k) = f(k)y'ra/ Ek- The equation 
obtained for f{k) has the relativistic form (^. It contains the factor 1/ek' in the integration volume and the kernel 
K-a' /r2 ik,k') — tj{k, k')^ek£k' / "m- For k' / k — j fixed the latter has the asymptotics (||) with the constant a' instead 
of the function An. For k' fixed, fc ^ cx) it decreases like 1/Vk, and analogously for k fixed, k' —* 00. 

The kernel K22, when /c — > 00 and k'/k = 7 fixed, has also the asymptotics (^, however with an unbounded 
function A22, in contrast to An (see below). When /c — > 00 and k' fixed (and vice versa) it tends to a constant. 
Therefore it always dominates over Kn and over K_^, i^'2{k,k'). 

To disentangle the two different sources of collapse, we first consider the one channel problem for the component 
/i with the kernel Kn, and we remove the second equation from (^. We analyse the domain of k'/k — 7 fixed, 
where Kn has its maximal asymptotics values. The 9,9' domain is finite and the function An{9,9',^) has no 
singularities in 9,9'. Therefore we majorate it by its maximal value. If the stronger, majorated kernel results in 
stable bound states, the exact one results in stability too. This method to analyse the cutoff dependence is equivalent 
to applying the sufficient condition of stability proposed in ||l^. The inspection of (|lO|) shows that for fixed 7, the 
maximum of An is achieved at 6* = 9' , where it has the form: An{9 — 9',j) = a'^/^ independent of 9. Note that 
^22 (^ = = 01! I ^ is unbounded when 7^0. Then we majorate the function Anif) = 0',l) by its maximal 

value relative to 7: An""^ = a', which is evidently achieved at 7 = 1. With this value of An, the kernel (||) exactly 
coincides with the kernel K_a'/r^{k,k'). As well known |Q, with the potential —a'/r'^, the binding energy does 
not depend on cutoff if a' < 1/(4to), what restricts the couphng constant to a < 7r/2. If a' > l/(4m), the system 
collapses, what is manifested by the fact that the binding energy tends to —00 when k^ax 00. In this system, there 
exists a critical value ac of the coupling constant, below which the binding energy is stable. Majorating the kernel, 
we underestimate Uc- We calculated also this value p^ , not majorating the function An = ct'\/l, but taking into 
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account its dependence on 7. In this way we find ac = tt, instead of 7r/2, when was replaced by 1. Because of 
majorating the kernel in the variables 0,9', this value should be still smaller than the true ac but it is in the range 
3 < Qfc < 4 we found numerically (see next section). 

In the two-channel problem, the kernel dominating in asymptotics is A'22- In the case J = it is positive and 
corresponds to repulsion. Because of that, this channel does not lead to any collapse. This repulsion cannot prevent 
from the collapse in the first channel (for enough large a), since due to coupling between two channels the singular 
potential in the channel 1 "pumps out" the wave function from the channel 2 into the channel 1. So, in the coupled 
equations system ^ the situation with the cutoff dependence is the same as for one channel. 

Let us now consider the state J — I, a = 0. The asymptotics of the kernel K22~^^ is the same than —K22~^\ it 
is negative and corresponds to attraction. Since it always dominates over K_a' /r^{k, k'), this attraction is stronger 
than in the —a'/r'^ potential. Therefore it results in a collapse for any value of the coupling constant. In the paper 
1^, this situation corresponds to the two-fermion state described by the components <I>^~,<I>^+. 



IV. NUMERICAL RESULTS 
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FIG. 1: Cutoff dependence of the binding energy in the 
J = or (l-|-,2— ) state, in the one-channel problem 
(/i) , for two fixed values of the coupling constant below 
and above the critical value. 



FIG. 2: Cutoff dependence of the binding energy, for 
J = (l-|-,2-) and J = l,a = (l-,2-|-) states, in 
the two-channel problem (a = 1.184). 



The preceding results are confirmed by numerical calculations. The constituent masses were taken equal to m=l 
and the mass of the exchanged scalar /x=0.25. 

We first present the results given by the single equation for fi with kernel Kn in the J — case. We have plotted 
in figure ^ the mass square of the two fermion system as a function of the cutoff k^ax for two fixed values of 
the coupling constant below and above the critical value In our calculations, the cutoff appears directly as the 
maximum value k^ax up to which the integrals in (^ are performed. One can see two dramatically different behaviors 
depending on the value of the coupling constant a. For a = 3, i.e. a < ac, the result is convergent. For a = 4, i.e. 
a > ac, the result is clearly divergent. decreases logarithmically as a function of k^ax and becomes even negative. 
We would like to notice that this divergence is not associated with the non decreasing behavior of the K22 kernel but 
with the existence of a critical value of the coupling constant separating two dynamical regimes. This property is due 
only to the large k behavior of Ku. Though the negative values of Af^ are physically meaningless, they are formally 
allowed by the equations (|l|) and (|^). The first degree of M does not enter neither in the equation nor in the kernel, 
and crosses zero without any singularity. The value of ac does not depend on the exchange mass 11. For n <^ m, 
e.g. /i w 0.25, its existence is not relevant in describing physical states since any solution with positive M^, stable 
relative to cutoff, corresponds to a < ac. For fi ^ m one can reach the critical a for positive, though small values of 
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We consider now the full Yukawa problem as given by the two coupled equations (0). In figure || are displayed the 
variations of for J = or (1+, 2—) and J = 1, a = or (1—, 2+) states as a function of the cutoff kjnax- The value 
of the coupling constant for both J is a = 1.184, the same that in Fig. 2 of below the critical value. Our numerical 
values are in agreement with the results for the cutoff A < 100 presented in this figure but our calculation at 
larger k,nax leads to different conclusion for the J = state. We first notice a qualitatively different behavior of the 
two states. In what concerns J = 0, the curve becomes flat when kmax increases, - with a 0.1% variation in 
when changing kmax from 50 to 600. We thus conclude to the stability of the state with J = 0, as expected from our 



analysis in sect. [II 



On the contrary, for J = 1, a = the value of M (kmax) decreases faster than logarithmically what indicates - as 

found in 1^ - a collapse. As mentioned above, the asymptotics of the i^22~^^ kernel is the same as the K^i~°^ one 
but with an opposite sign, i.e. it is attractive, what leads to unstability for any value of a. We found the same result 
when solving the J = equations with the opposite sign of • 



V. CONCLUSION 



The Light-Front solutions of the two fermion system interacting via a scalar exchange have been obtained. We have 
found that the J = - or (1+, 2—) - state is stable (i.e. convergent relative to the cutoff kmax — > 00) for coupling 
constant below some critical value, in a way similar to what is known in non relativistic quantum mechanics for the 
—a'/r^ potential. In this point, our conclusion differs from the one settled in where it was stated that the integrals 
in eqs. (|l|) diverge logarithmically with cutoff. Above the critical value the system collapses, what manifests as an 
unbounded cutoff dependence of Af^. In the J = l,a = - or (1— ,2+) - state the system is found to be always 
unstable, in agreement with 

The origin of this unstability for the J = state differs from J — I. We have found that the iir22~"' dominating 
kernel does not generate a collapse because it is repulsive. The unstability in this case is related to Kn. 

These results should be taken into account when carrying out the renormalization procedure. The explicitly 
covariant LFD may be efficient for solving this problem, like it has proved to be fruitful for analyzing the Yukawa 
model. 
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